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Cosmology in extended theories ol gravity is considered assuming the Palatini variational principle, 
for which the metric and connection are independent variables. The field equations are derived to 
hnear order in perturbations about the homogeneous and isotropic but possibly spatially curved 
background. The results are presented in a unified form applicable to a broad class of gravity 
theories allowing arbitrary scalar-tensor couplings and nonlinear dependence on the Ricci scalar in 
QQ ■ the gravitational action. The gauge-ready formalism exploited here makes it possible to obtain the 

' equations immediately in any of the commonly used gauges. Of the three type of perturbations, the 

, main attention is on the scalar modes responsible for the cosmic large-scale structure. Evolution 

^SJ ■ equations are derived for perturbations in a late universe filled with cold dark matter and accelerated 

by curvature corrections. Such corrections are found to induce effective pressure gradients which are 
problematical in the formation of large-scale structure. This is demonstrated by analytic solutions 
. in a particular case. A physical equivalence between scalar-tensor theories in metric and in Palatini 

formalisms is pointed out. 

0^ ■ 
^ ■ 

> ■ I. INTRODUCTION 

(N : 

, The observed acceleration of the universe [l|, 0, 0| might indicate that our understanding of gravity breaks dovifn 
■ at cosmological distances. In the context of general relativity, what is needed is some sort of fluid with large negative 
' pressure to accelerate the expansion. In its simplest form this fluid is the cosmological constant 5] , equivalent to a time- 
I independent vacuum energy density. But there is a problem with the cosmological constant, namely that according to 
I present knowledge of quantum field theory, one would expect this vacuum energy density to be larger than observed 
. by an enormous factor. One approach to this problem is to assume that due to a symmetry the cosmological constant 
CIh' is zero, and in its stead a dynamical fluid, generically called dark energy 0|, provides the negative pressure. A plethora 
Q , of different dark energy models have been proposed, but none without some problems of its own. 
^ ■ As an alternative to energetics of unknown fluids, it has been proposed that the cosmic speed-up could stem from 
I modifications to general relativity 0, H, H, 113, IHi- Such modifications can come into play when the gravitational 
d • action contains other curvature invariants apart from the standard Einstein-Hilbert term. Then general relativity 
can be considered as a limit of a hypothetical more general theory. In fact suggestions for such a theory can be 
found from fundamental physics. Quantization on curved spacetimes has been found long ago to require extension 
of the Einstein-Hilbert scheme by addition of higher-order curvature termsfT2l|, and it has been also shown that the 
dominance of these terms at high energies could be the cause of the early infiationary period in the universe[13, 14, 15j]. 
Similarly, corrective terms that become iniportant at small curvature can lead to late-time effects seen in the cosmic 
expansion as an effective dark energy [H, [l^, [H, [H, HO, HU, HI]- It has been shown that such curvature terms can 
arise from compactification of time-dependent extra dimensions in string /M-theorv [23| . 

Once the gravitational action is nonlinear in R, the question which variational principle to apply becomes relevant. 
The Palatini variation of a nonlinear gravity action leads to a different theory than the metric variation. The metric 
variation of extended gravity theories result in fourth order differential equations which are difficult to analyze in 
practice. The Palatini formulation, in which the connection is treated as an independent variablef23| is more tractable 
than the metric one, and it can also in general exhibit better stability properties (25i, i26i]. since it yields the modified 
field equations as a second-order differential system. Mathematical convenience does not of course prove that the 
Palatini variation would be the fundamentally correct procedure. However, this possibility might be interesting also 
according to some theoretical prejudices. The second-order nature of the Palatini formulation is conceptually more 
reconcilable with better-known physics than the metric alternative, where the action in the beginning contains second 
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derivatives of the metric, and in the end one has to specify initial values up to third derivatives to predict the evolution 
of the system. The doubling of the variational degrees of freedom in the Palatini formulation has an analogy with 
the Hamiltonian mechanics, where the coordinates and momenta of particles are treated as independent variables [27j. 
On the more speculative side, it is also interesting that the Palatini scheme of gravity can be recovered in unification 
of general relativity with topological quantum field theory .28.]. 

Nonlinear Palatini gravity may also pas s the Solar system tests [29|. While this result has been questioned [sol Isii] . 
the most recent discussions support it|32l. [ssl . [3^ . Cosmological background solutions have been given for several 
choices of the gravitational Lagrangian in the Palatini formalism, including (with some constants a and /3) / (R) = 
-R^dl, f{R) = alog(/3i?}jp, f{R) ^ R + aR-l^^ [37|, f{R) = R + aR'^^, f{R) = R + a/R + I3R^^ and 
f{R) = R + a/ R + /3i?^[3J]. Such solutions have been tested against observational data, confirming that the late 
accelerating expansion history can be produced viably in these models [1^ H^, |4l|. This suggests that one should 
investigate the cosmological implications of these theories beyond the background order. While a wide variety of 
models is able to produce a scale factor evolution consistent with observations, considerations of the cosmic microwave 
background anisotropies and large scale inhomogeneities in cosmic structure may distinguish between different physical 
assumptions underlying similar or even identical expansion histories [4^. |43L H. l45l| . Recently it has been shown that 
early inflation might also be modelled by nonlinear gravity in the Palatini approach 3J| , and it will be interesting to 
see the primordial perturbation spectra such models predict. 

The aim of this paper is to set up perturbation equations for generalized gravity in the Palatini approach. Cosmo- 
logical perturbation theory 46, 47, 48] of modified gravities in the metric approach has been rigorously formulated in 
the gauge-ready form by Hwang and Noh[4^ [s^, [sij. Recent extensions of these formulations encompass also kinetic 
theory [53], tachyon condensation [E^ and string corrections [s^. The pioneering works found in these references pro- 
vide an arsenal of efficient methods on which we draw to derive and analyze the Palatini versions of the perturbed field 
equations. In particular, we will exploit the gauge-ready formalism where the equations are instantaneously obtainable 
in any of the convenient gauges employed for cosmological perturbations in the literature. Since in the general case 
it is not clear a priori which gauge choice is the most suitable to simplify the analysis, leaving the temporal gauge 
unfixed in the equations gives them very advantageous adaptability to different situations. 

In Section II we will review briefly the Palatini variational principle. We also point out an equivalence between 
metric and Palatini formulations of scalar-tensor theories, which is verified in detail in Appendix A. In Section III 
we identify the degrees of freedom in the cosmological metric and the energy- momentum tensor, and derive the 
corresponding field equations for the three types of perturbations, scalar, vector, and tensor. In Section IV we 
discuss structure formation in f{R) cosmologies, in particular curvature corrections that become important in the 
late matter-dominated universe. Section V contains discussion of our results. 



II. PALATINI APPROACH TO GENERALIZED GRAVITY 



We consider gravity theories represented by the action 



d 



xy'-g 



(1) 



Here ... are some matter fields. In the Palatini approach one lets the torsionless connection f vary independently 
of the metric. The Ricci tensor is then defined solely by this connection. 



whereas the scalar curvature is given by 
The matter energy-momentum tensor is defined as 



R = g^'R^, . 



2 &{^C„ 



(2) 
(3) 

(4) 



Hereafter we drop the sub- and superscripts m, which stand collectively for all matter except the scalar field (f>. We 
consider the second and the third term in the action ([1]) as the Lagrangian for 0. Then the scalar field energy- 
momentum tensor 



(5) 
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is found by a similar variation as Eq.(|4]). The field equations which follow from extremization of the action, Eq.(IT]), 
with respect to metric variations, can be written as 

FR^-^fSi^ = Ti: + Ti:(^\ (6) 

where we have defined 

F = df/dR. (7) 

In general relativity, / = R/8ttG or {R - 2A)/8ttG, and F = I/SttG. 

Note that the covariant derivative V in Eq. ([5]) is taken using the Levi-Civita connection (i.e., the Christoffel 
symbol) of .g^^,, 

r)37 = 2 (5A/3,7 + 5A7,/3 - 5/37, ■ (8) 

Thus we have two different connections, the "Palatini" connection E^^, appearing as an independent variable in the 
action ([1]), and the Levi-Civita connection E^^ derived from the metric g^^,, and two different covariant derivatives, V 
and V, corresponding to these two connections. We denote by Rfj.,y{g) = E^^^ — E"^ ,^ + E";i^E^^ — T'^x^au the Ricci 
tensor constructed from the Levi-Civita connection to distinguish it from the Ricci tensor i?^^ of Eq. Likewise, we 
write R{g) = g^'^ Rfiuig) to distinguish this scalar curvature, derived solely from the metric, from the scalar curvature 
of Eq.([3|) appearing in the action 

By varying the action with respect to E^^, one gets the condition 

[V^g^-'F] = 0, (9) 

implying that this connection is compatible with the conformal metric 

ij,,, EE (10) 

This connection governs how the tensor R^^ appearing in the action settles itself in order to minimize the action, but 
it turns out that the connection of Eq.® determines the geodesies that freely falling particles follow, since energy 
momentum is conserved according to this connection [s^, 

V^T,^ = 0, (11) 

whereas in general V^T,^ ^ 0. Therefore one may interpret the Levi-Civita connection as the gravitational field as 
in Ref.[56]. In any case it is convenient to write the modified field equations in terms of the connection E. One finds 
that the Ricci tensor is 

R,. = R,A9) + - ;^(VmV.f) - ^-L_}-g^,uF. (12) 



The curvature scalar and Einstein tensor follow straightforwardly. 



R = R{g) 



[n - 2)F 



-2UF +—{dFf 
F 



(13) 



G,. = + |!i-lli^(V,J^)(V.^^) - 1 (V,V. - g,,U) F ^^-^^l-^g^,[dFf . 



(14) 



Erom now on we set the spacetime dimension to 7i = 4 and use units SttG = c = 1. One useful way to write the field 
equations is in the form of Einstein gravity plus corrections: 



Gi{g) 



(1 ' F)R^,{g) 



A(vMF)(V.F) + (V''V.F) + i 



{f-R) + ii-j)^F + ^idFr 



St 



(15) 



Due to Eq. (|lip and the Bianchi identity, the correction terms are covariantly conserved according to E. In the metric 
formalism, analogous corrections appear even in vacuum, and it is natural to treat them as an effective additional 
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fluidfiQ*]. However, here these corrections can be expressed as functions of the trace of the matter energy-momentum 
tensor T = g^'^T^^, and the scalar field and its derivatives, and one can view Eg.p^ as general relativity with 
nonstandard matter couplings. Thus the whole right-hand side may be regarded as an effective matter energy- 
momentum tensor. In vacuum it reduces to a cosmological constant[57], and in the case of conformal matter (i.e. 
T = [1^) to the usual T^^. One finds R in terms of matter from the trace of the field equation, Eq.®, 

FR~2f ^uj{d(l)f -W + T. (16) 



We will refer to this central relation as the structural equation [59|, |60]. 

Now consider some specific cases. Perhaps the simplest nontrivial example is a pure scalar-tensor theory, for which 
/ = F[(p)R. Then the structural equation becomes 

R^-"-l^^. (17, 

We remark here that this theory is physically equivalent to a similar one in the metric formulation with a modified 
kinetic term of the scalar field. Namely, with arbitrary F{(j)) and Lo{(j))^ the kinetic coefficient ^metric for the scalar 
field in the corresponding metric formulation of scalar-tensor theory is 

u;„,,Hc(0)=a;(</>)-^^. (18) 

We show in Appendix A that the metric and Palatini approaches yield exactly the same field equations and equations 
of motion when this rescaling of oj is taken into account. Therefore no qualitatively new features, regardless of the 
scalar field couplings, are introduced in the Palatini approach whenever the gravitational action is linear in R. This 
correspondence, Eg. p^ . enables one to extend the results of studies of scalar-tensor gravity^] within the metric 
formalism to the Palatini framework, and may motivate some previously less-investigated forms of these theories. 

Consider the case without a non-minimally coupled scalar field, / = f{R). We assume that there exists a solution 
R = R{T) to the structural equation FR — 2/ = T. For example, if f — R — ao/{5R), then 

RiT) = -i (t+ v/r2+4Qo) . (19) 

If f{R) ^ R + ctoR^, we have R{T) = — T as in general relativity. If / = R^, the theory admits only vacuum or 
conformal matter. 



III. COSMOLOGICAL EQUATIONS 
A. Definitions 



The line-element in the perturbed Friedmann-Robertson- Walker (FRW) spacetime can be written as (see e.g. 
Ref. [i^l for an only slightly different notation) 



(r,) {- (1 + 2a) dif - 2 (/?,, + h) di^dx' -f [gf ^ + 2 (fflf <^ + 1\^J + C(,|j) + /ly)] dx'dx^^ 



(20) 



We characterize the scalar perturbations in a general gauge by the four variables a,P,ip,j. Vector perturbations 
introduce four more degrees of freedom, the divergenceless 3- vector fields bi and c;. Gravitational waves are described 
by the two free components of the symmetric, transverse and traceless 3-tensor hij. We have thus decomposed 
the ten independent components of the symmetric tensor Sg^i, into three types of perturbations according to their 

transformation properties under spatial rotations. The comoving spatial background metric g^j reduces to 5ij in a 

(3) 

flat universe. The vertical bar indicates a covariant derivative based on the Levi-Civita connection of g^^ . This metric 
is used to lower and raise spatial indices i,j,k... of the perturbation variables. 
The components of the energy-momentum tensor for an imperfect fluid are 



= -ip + dp), T^ = -{p + p) (y,, + v^''^) , = (P + Sp)S) + 



(21) 
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Here p and p are energy density and pressure, and w, w*^"^ are the scalar and vector velocity perturbations"'^, respectively. 
Background quantities are denoted with an overbar, which we will usually omit when unnecessary. The isotropy of 
background does not allow anisotropic stress except as a perturbation. This we decompose into the scalar, vector and 
tensor contributions as 

n,, = ( n[^^ + -^n^A + n|''), + n|*^ , (22) 

*J I |ij ' 3 / ' (i\3) «J ' ^ ' 

where A stands for the three-space Laplacian based on the Levi-Civita connection of g^J^ . The vector n^"-* is diver- 

genceless and the tensor lif^ is symmetric, transverse, and traceless. The four scalar degrees of freedom for the fluid 

perturbation are therefore (5, 5p^ v, and II^''-', independent components of the divergenceless vectors w^-^'' and n[^'' sum 

up to four and the tensor describing gravitational waves II^*'' has two independent components. 

Some of these degrees of freedom are due to arbitrariness in separating the background from the perturbations. In the 
gauge-ready formalism [50j one deals with these gauge degrees of freedom by noting that the homogeneity and isotropy 
of the background space implies invariance of all physical quantities under purely spatial gauge transformations. 
Therefore one can trade (3 and 7 to the shear perturbation 

X = a{(3 + i). (23) 

where an overdot means derivative with respect to the conformal time ?/. Since both f3 and 7 vary under spatial gauge 
transformation, they appear only through the spatially invariant linear combination x in relevant equations. In 
addition, one can define the perturbed expansion scalar 

n=-iHa-ip)-^X, (24) 

Here H is the Hubble parameter defined with respect to conformal time, i.c, the usual Hubble parameter multiplied 
by the scale factor a. The variable k is a convenient linear combination, the use of which simplifies some equations, 
but it is not linearly independent of other perturbations. Only three of the variables o;, (£, x f^^id k are independent. 
All of them can be physically interpreted as perturbations of the normal- frame vector [50j. The advantage of using 
this set of variables is based on the fact that they are spatially gauge- invariant. Thus, writing equations in terms 
of them, one can conveniently fix the temporal gauge by just setting one of these metric perturbations to zero. The 
synchronous gauge, corresponding to a = 0, is an exception where the gauge mode is removed only up to a constant. 
For example, the longitudinal (also called the Newtonian or the zero-shear) gauge is the one where x — 0- One 
popular gauge is the comoving one, where, instead of setting any of the metric perturbations to zero, one sets the 
fluid velocity perturbation v in Eq. (PTjl to v — . Suitable linear combinations of the above gauge conditions can be 
considered also. 

Similarly, we will exploit the spatially gauge-invariant variable 

= h + c, (25) 
to characterize vector perturbations of the metric. The tensor perturbation hij is gauge-invariant by construction. 



B. Background 



The modified Friedmann equation corresponding to Eq.® is 



3H^ = - 



9 ^ )9 



a^V - 3HF - —F^ - — (/ - FR) 
AF 2 ' 



- 3K, 



(26) 



where K is the curvature of the background space. One should be equipped with a solution to the structural equation 
Eg. pop , with which to replace the scalar curvature R appearing both explicitly and implicitly (through /, F and 



^ Note that we use the covariant velocity perturbations [6^ , sometimes denoted as d = a{V — (3), and i)'"' = a(V^'"' — bi). 
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their time derivatives) in Eq. p6|) . For a review of background solutions in several f{R) cases, see[63|- An implicit 
expression for R follows from Eq. (fT51) : 

a^R = 6 (^H + + k'^ + J (^F + 2HF - ^F^^ ■ (27) 
The matter continuity equation is, as usual, 

p + 3H{p + p) = 0. (28) 

Individual matter species, if minimally coupled to gravity and without interactions with other matter species, also 
satisfy this equation. The Klein-Gordon equation for the scalar field is 

4> + 2H^+^^' + ^{2V-f)^^^Q. (29) 



C. Field equations for scalar perturbations 



First we present the equations that govern the evolution of scalar perturbations defined in the subsection IIIA. 
From now on we consider variables in the Fourier space. The transformation is simple since at linear order each 
/e-mode evolves independently. 

The energy constraint (Gq component of the field equation) in generalized gravity is 



2H 



F 
F 



+ {QK-2k^)ip + 



^ujcf,^ + — +ZHF\ a 



i J - aHp - w^5^ - i + a" {2V - /)^, 



3 (H^ + 



4F2 2^^^ 



SF 



3£ 
2F 



3H]SF 



(30) 



and the momentum constraint (G^ component) is 



,,n „ . 1 3F 3 



{p + p) - + uj 



H 



3F 
2F 



SF + 5F 



(31) 



The shear propagation equation (G* — j'^jG'^ component) reads 



H 



- - 

F / a 



1 (a^nf^) + 5F^ . 



(32) 



The Raychaudhuri equation (G^ — Gg component) is now given by 



2ak 



AH 



F 
F 



6 H -H' 



F_F^ 
F^ F^ 



F 4 



2k^ 



F 

a + 3— d 
r 



f' 



{5p + 3(5p) + 4w0(50 + 2w>^ + a (/ - 2V") 



F 



F ■ 



+ 3— r - a^R + QH + \ 5F - 6—SF + 36F 



■p2 



F 



(33) 



The scalar curvature perturbation can be found by combining previous equations or directly by using Eq. (ll3|) . 



a^5R 



+ 



-2ak- I 87J + 3|J j aK + 3 



F2 F F 



F 



F^ F 



2\H^-H--\+ — - —H + -k^ 



6F+ \2H 



F 

f' 

F 



F2 F 
5F + 5F 



a + 2 (2fc2 -QK)<f- 3^a 

r 



(34) 
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Also redundant with Eqs. p0l - [33|) is the perturbed version of the structural equation (G" component) 

R5F + FSR - 26 f = ^uj(j) (s^ ~ + 2(/ - 2V),^S(I) -Sp + 3Sp. (35) 



We have Sf = FSR + f.^Scf), and similarly SF = FuSR + F^^S4). However, these functions can also be expressed 
in a form proportional to matter perturbations using the structural equation Eq. (|16p . Consider the earlier example 
f = R- ao/iSR), Eq.lUni). There we have 

ao \ f 1 T 

4q:o 



SF 



{-Sp + 3Sp) 



(37) 



It is straightforward to obtain the derivatives SF and SF from Eg. ([57)1 . In this manner SR can be related to ST, 
and one can consider the right-hand sides of the perturbed field equations as matter sources in the modified gravity. 
However, one possibly appealing gauge condition is to set SF = 0. This eliminates the complicated effective matter 
source terms in the right hand side, while leaving some /-dependent terms to modify the evolution of the metric 
perturbations. In the / = ,f{R} case this is the gauge where ST = 0. 

The equation of motion for the scalar field is the same as in the metric formulation, 









+ 2H^S(j) + 







II 



Ijj 
2lo 



(f) (a + an) 



LU 



(0 + 2i7) 



2uj 



2(j) + H(t) + 2—4)' 

UJ 



(38) 



However, now the interpretation of R appearing on both sides of Eq. 
matter species obey the standard continuity and Euler equations. 



is different. The minimally coupled individual 



Sp ^ 3H {Sp + Sp) ^ (p + p) {-kv + aK- iHa) , 



(39) 



1 2 fc^ — ZK 
— [a^{p+p)v\ ^k{p + p)a + k5p- H(^), 

since matter fields in the action, Eq.ljT]), decouple from the connection. 



(40) 



D. Vector and tensor perturbations 

Vector and tensor modes are simpler in an FRW spacetime than scalar perturbations. This is true especially in 
f{R, (p) gravity, since neither R nor (j) has vector or tensor components, regardless of whether we are in the Palatini 
or metric formulation. Therefore we can employ the analogous results found in the metric case, see for example 52]. 
The matter quantities in this subsection do not involve contributions from the scalar field, since it cannot generate 
vector or tensor perturbations. 

The equations governing the evolution of rotational perturbations are 

k"^ -2K \ , ^ M 

-—^^: = -{p + p)vr, (41) 



a^{p + p)v'f^ 



fc2 - 2K^ 
2fc2 



-H 



{v) 



(42) 



The first is the field equation (the component), the second one the conservation equation. The only difference to 
general relativity here is the appearance of the prefactor 1/ F modulating the amplitude of 'Si. For example, one can 

see that the angular momentum of a perfect fluid ^ a^(p -f p)f|"'' is conserved in an expanding universe. 
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The gravitational wave is governed by the equation 

/i; + (^2H + 1^ + (fc2 + 2K) h) = ^Tlf. (43) 

Modifications to general relativity appear here as an additional damping term and a similar modulating of the source 
term as for the rotation. In the absence of anisotropic stresses, only the former has an effect. For a flat universe, an 
integral solution exists for the superhorizon scales, 

h^(r^,k)^A^^B)f -L^dr,, (44) 

where Aj and Bj are constants for each fc-mode. Thus only the decaying solution is modified when / ^ R. 



IV. STRUCTURE FORMATION IN f{R) COSMOLOGIES 
A. Inhomogeneities in late universe 

Now we are in a position to consider formation of structure in cosmologies derived from the Palatini approach 
to generalized gravity. As discussed in the introduction, it is of particular interest to uncover whether curvature 
corrections driving the late acceleration will predict testable features in the CMB or large scale structure. Therefore 
we will consider f{R) theories in a flat universe dominated by pressureless cold dark matter. There is no dark energy, 
since we consider modifled gravity as its alternative. 

The evolution of matter perturbations is governed by the equations ([39]) and (j40|) . which for pressureless matter 
reduce to 

S = —kv + aK — 3Ha, v = —Hv + ka. (45) 

We will work first in the uniform density gauge (indicated by the subscript s), where 6s — 0. Now also Sfs and 
SFs together with their derivatives vanish, considerably simplifying the analysis. (This of course does not mean that 
the matter perturbation disappears, it is just carried by other perturbation quantities, as the physics is completely 
independent of gauge choice. We later transfer to the comoving gauge, where the interpretation of the results is more 
straightforward.) 

In the uniform density gauge, manipulation of the field equations (jSOHSSp with the help of the conservation equations 
(|45p then yields an evolution equation for the non- vanishing matter velocity perturbation vg, 



1 



F{2FH + F) 



-2F^(H^ + 2H) + 2F^ - FFH - 2FF 



i.2 

GF^HH - 2F^H + FF{H + — ) + 2FFH 

3 



Vs 



(46) 

By solving this relatively simple (depending on the form of f{R)) differential equation, one can easily find also the 
metric perturbations, since they are related to vs and its derivatives. These solutions can also be related to solutions 
in any other gauge by using the gauge transformation properties of the relevant variables. In fact, although the 
evolution equation (|46p is particularly tractable in the uniform-density gauge, it is difficult to physically interpret the 
significance of modifications to general relativity from this equation. 

Therefore we transform this result into the more intuitive comoving gauge (subscript „), which is defined by 
Vy — 0. We cannot of course find a nontrivial equation for v^, but will instead consider the comoving gauge density 
perturbation 5y. By looking at Eg. ([45]) one sees that this gauge also coincides with the synchronous one at the 
late matter-dominated stage of the universe, regardless of the form of f{R)- By using their gauge transformation 
properties, we find a convenient relation between the perturbation variables in the different gauges, 

(5„ = 3Hvs/k. (47) 

Using this we can recast Ea. (|46| into the form 

S,, = ^ \ -3H \2FH(FH^ + F) - 2F^H + FF(-2H + H^)] 5^ 

3FH^{2FH + F)^ I ^ ^ 

QF^H^iii - 2HH) + 6F^H{H^ - H) + FF{3HH - 6H^ - H^k^) + 6FFH{H - H^) (5„|. (48) 
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We have checked that the form of this equation does not change when the assumption if = is relaxed. 
For the case of general relativity with a cosmological constant, f{R) — R — 2 A, Eq. (H51) reduces to 




S, = -HSy+\—-2H]S,. (49) 



If the cosmological constant and cold dark matter arc thought of as components of a two-component cosmic fluid, 
the equation of state of this total fluid is wt = —(1 + a^^Po/A)^^. The density contrast (5j can then be written as 
(5t = (1 + wt)Sv. Inserting these in Eg. p^ gives the familiar 

Sj. = {6wT - 1)H5t + ^[1 + 8wt - ?>wI]H'^5t. (50) 

Notice that this equation is scale-independent, whereas Eq. (|48p involves a fc^ -term, which stays small only at very 
large scales. 

Consider first this large scale limit of Ea. (|^5)) . There both the "friction" and the "source" term in Eg. ([50)1 are non- 
trivially modified. Similar effects appear when the cosmological constant is replaced by dynamical dark energy [3 [g^. 
Here these modifications depend on the specific form of f{R), and one might recover an evolution of Sy leading to a 
large-scale matter power spectrum consistent with observations, by choosing a suitable function of R. 

However, there could still be difficulties with the small-scale structure in these models. This is because the gradient 
term in Eq. (|48p acts as the dominant source inside the horizon. The gradient appears since the curvature corrections 
induce effective pressure fluctuations in the inherently cold dark matter. This may have been anticipated from the field 
equation (jisp , where it is seen that the modified Palatini gravity is coupled to derivatives of matter energy momentum. 
Since the cosmological background is homogeneous, modifications to the Friedmann equation stem only from the time 
variation of the background density of matter. However, the evolution of inhomogeneities in the universe is inevitably 
affected by the response of the modified Palatini gravity to spatial variations in the distribution of matter. This is 
why we have found here an effective sound speed for cold dark matter. 

It might be useful to compare the situation with dark energy models. A cosmological constant does not fluctuate, 
therefore the absence of a gradient term in Eg. ljSOp . On the other hand, quintessence fields fiuctuate but their gradients 
are negligible, since those fields are very smooth inside horizon, although this may not be the case if the quintessence 
field is coupled with dark matter [5^. The models unifying dark energy and dark matter into a single component, 
are plagued by a nonzero sound speed at small scales [43, [6J|. These unified models are somewhat analogous to the 
case studied here, since we assume that due to the properties of non-standard gravity, the accelerated expansion of 
the universe, usually ex plai ned by dark energy, is in fact driven only by dark matter. Stability of modified gravities 
in the metric approach 6^ 6^ has been explored in de Sitter space using covariant and gauge invariant analysis of 



cosmological perturbations [67|, |68|, but it remains to be studied whether similar effects in the evolution of matter 
inhomogeneities as seen here will appear in the metric approach. 

Here, at sufficiently large fc, the small scale evolution of is determined by 

.. F'{T)p 

2F(T) + 3F'(r)p^ ^^^> 

where we have expressed the gradient in Eq. ([48]) in terms of the background solution for F one finds from the structural 
equation. It is interesting to note that the gradient (or the effective sound speed) does identically vanish only in the 
case of Einstein gravity with a cosmological constant. For example, when f{R) = R — ao/{3R) one has 



Sv = - 



2x 



k^Sy, (52) 



_3(1 + 6a;) -{3 + 8a;) VI + 4a;. 

where the square bracket term is of the order one when x = aop~^ is, which holds when the curvature and matter 
energies are both relevant. Indeed the effective pressure fluctuation seems to become troublesome inside the horizon. 



B. A particular example: f{R) ~ R" 

We demonstrate the general considerations of the previous subsection with a specific choice for the nonlinear 
Lagrangian, /(i?) ~ R"', where n ^ 0,2,3^. Such a case has been studied previously [3^ [6^ and shown to predict a 



^ The limit n = 3/2 appears singular in our equations. This is because at that limit the deceleration parameter vanishes and the conformal 
Hubble parameter cannot be expressed in the form ~ l/rj. When n < 3/2, is positive, and when n > 3/2, is growing from negative 
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plausible background evolution at late times [4l|. In fact the background is simply described by a constant effective 
equation of state in this model. The Hubble parameter scales as ~ a^~'^/". It is then easy to write it with its 
derivatives in terms of conformal time, 



1 



(3-2n)2 7y2^ 



(53) 



This solution applies also in the metric formulation [70|. In addition, Palatini variation of Ricci squared gravity, 
/ ~ {W^ R^iv)^/'^ , gives the same expansion ratefooj. 

Here the scalar curvature is i? = 3(3 — n) H"^ / {2na^) . Since Eq. (|48| involves only F/F and F/F, we need only to 
know that 



F 



3(1 -n) 



HF, F 



3(1 -n) 
2n2 



(3 - 4n)H'^F. 



Substituting these in the general Ea. (|i5)) yields 

3. 



{n- -)S^ 



V 



3 (n-l){n 

h ^ — 

,2 ^ 



2) ■ 



n 



(54) 



(55) 



Inserting a power law ansatz S-^ ^ a™ one finds that in the large-scale limit fc = the solutions are m — — l/(6n) and 
m = 3/n — 2. When n is positive, the former is the decaying solution, and the latter gives the growing mode, which 
reduces to the standard matter-dominated solution 6 ^ a when n = 1. The growth rate reduces when n increases, 
reflecting the fact that the background expansion is sped up. In the limit n — 3/2, when a — 0, then also 5i, — 0. For 
larger values of n the large-scale inhomogeneities in fact begin to smooth out, instead of just their growth slowing 
down. 

In the small-scale limit k^r/^ 3> 1, and the gradient term drives the perturbations to exponential growth 



6y ~ exp ± 




(56) 



when 1 > n > 3, or to oscillate. 



exp ±i 



-krj 



(57) 



when n > 3. Since such features are not observed in the matter power spectrum, these models do not appear to be 
viable alternatives to dark energy. A quantitative study using existing observational data is under progress J Ij. 



V. DISCUSSION 



In the metric formulation nonlinear gravity theories result in rather untractable fourth-order differential equations. 
In the Palatini formulation the structure of the theory becomes interestingly different. Since there is an algebraic 
relation between the curvature scalar and the trace of the matter energy-momentum tensor, the Palatini variation 
yields second-order field equations. In Sec. I we discussed observational and theoretical motivations for studying these 
theories. In this paper we have demonstrated that even the theory of cosmological perturbations is solvable for the 
Palatini formulation of modified gravities. 

In Sec. II we presented the equations for linear cosmological perturbations in a general /(-R, 0)-gravity. We arranged 
them in a gauge-ready form, from which one can easily adopt the most suitable gauge for the particular problem at 
hand. We found that evolution of vector and tensor perturbations are modestly modified from general relativity, 
featuring some /-dependent prefactors. Field equations for scalar modes of perturbations are modified in a more 
intricate manner. In particular, nonlinear Palatini gravity features nonstandard matter couplings, originating from 



values towards zero. Our results are still well defined at the limit n = 3/2. On the other hand, n = 2 is a physically singular limit of 
the theory. This is not seen explicitly in our equations, since factors like l/(n — 2) belong to an otherwise irrelevant proportionality 
constant in the solution for / which we have rescaled away. We have also a singularity at n=3 [SSll . 
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perturbations in the trace of the fluid energy momentum tensor and its derivatives. Only in a specific gauge, 6T = 0, 
these couphngs ostensibly disappear. These equations can be almost trivially generalized to f{R, (p, X) gravity, where 
X = {d(j>)'^. For simplicity such a case was not included here, but curvature cou pling s of scalar field derivatives have 
been introduced as an approach to the cosmological constant problem fiol [60l \7% l73l . and it might be interesting 
to investigate implications of also these models to cosmological perturbations. 

In Sec. Ill we derived also equations governing the growth of inhomogeneities in the late universe undergoing an 
acceleration of expansion driven by nonlinear /(i?)-gravity. We found that an effective pressure gradient term appears 
in the evolution equations for cold dark matter inhomogeneities. This seems to be problematic for structure formation 
in such models. Analytical solutions were given for a simple example, f{R) ^ i?", but the problem seems generic 
to these nonhnear gravities. In the f{R) — R — ao/{3R) case, Eg. ([5^ imphes that the gradient term will drive the 
perturbations at sufficiently small scales. It might be difficult to tune the form of f{R) (unless it is trivially close to 
f{R) — R~ 2A) in such a way that the gradient term Eq. ([5T|) would be small enough inside the horizon not to affect 
the linear matter power spectrum. However, an extensive study of the constraints on the form and parameters of the 
function / in the action H]) is left for forthcoming studies. 



APPENDIX A: THE EQUIVALENT METRIC FORMULATION OF A SCALAR-TENSOR THEORY 

In n dimensions Eq. (fT8|) generalizes to 



i^metric{4') = w(0) 



(n-2) F{(p) 



Let us therefore consider the action 



S ■ 



d''x^{ -F{^)R{9) 



21 



in -2) F{<j)) 



(90)2 -y(0)+A„(.g^„ $,...) 



Variations of the metric lead to the field equations 

FR'i{g) - \fR{9) = n + (V^V, - <5^;n)F 



UJ - 



{n - 1) F'"^ 
(n-2) F 



{V^^<i>){VA)-5>i 



[n - 1) F'^ 
(n - 2) ~F 



which can be rewritten as 
1 



(Al) 



(A2) 



{dcj)f + V 



(A3) 



F" 



[n - 1) F'"^ 
{n-2)~F 



(V^<^)(V,(/.) + f^'(V''V.0) + <5^: 



-F" 



(n - 1) 

2(n-2) F 



{d(j)f - F'Uq 



,(A4) 



where T^'^'^'^ is given by Eq.®. 

Now let us vary the action of Eq.ljT]) when / = F{(j))R. We then obtain the Ricci tensor from Eq.([T 



R>i = Rig)^ + ^ 



{n-l)F'' 



in-2) F 
and the Ricci scalar from Eq. (fT3|) 



F" 



(V^0)(V.0) + F'(V^V.0) - J^ 



[F"id(bf +F'a^ 



(A5) 



R = Rig) 



(n-l) 1 
(n -2)F 



pl2 

— - 2F" ] {d(l)f - 2F'aq 



(A6) 



It is then straightforward to plug these into the field equation ([6]) . The result is exactly Eq. (jA4[) , which verifies that 
the scalar-tensor theory F{(f>)R with the Palatini variation is equivalent to the same theory with the metric variation 
and the kinetic term rescaled according to Ea. (|Al[) . The equations of motion for the scalar field are built into the 
field equations [5^, they also coincide in these equivalent theories. 
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APPENDIX B: SCALAR PERTURBATIONS OF THE METRIC 



In this Appendix we present some formulae which have been used in the derivation of the field equations. Here 
all covariant derivatives and curvature variables correspond to the Levi-Civita connection of the metric g of Ea. (j20p . 
These results do not depend on the variational principle which one uses. For a more extensive collection of useful 
formulae for metric perturbations, see appendices in Ref . [75| . 

The Levi-Civita connection in the perturbed FRW spacetime, Eg. ([^0)1 . is 

r[!o = + a , r", = (a - , ro, = g^f [Hil - 2a + 2^) + ^] + {2Hj + 7 + , 



J3) 



3k 



Useful contractions of these are 



1J.0 



fii ki 



4iJ + d + 393 + 7' ^ 
Covariant derivatives of a field S,{x,t) = ^{t) + (5^(5;, t) are then 



fc(3) 



a + + 7l\. 



1 

a 



5]{2Ha-^)~-x'' 

-I - 2hI -5\- 2H5^ + (J^l'-' fc + (2I + HCja + |(d + ok) . 
The Ricci tensor can be calculated using the "hatless" version of Eq. ^ , 



with Eqs. (|BT|1 . 

a'^(.9)j, 



3H~ak~ 2HaK + 3{H^ - H)a - Aa , 



= 2 



Ha-(p + {H - H^)(3 



2[-Ha + a]^^-R{g)ll^P'\ , 



6)iH + 2H^)+R{g)f^ +S) 



1 H 

-a-ip+-x-\ X 



(p~H{a- 5ip) - 2{H + 2H')a - 2i? + A(-^ + — x) 

a 



H 



The curvature scalar is 

a^Rig) = 6{H + H'^) 



-2 ak + AHan + 3{H - H'^)a + + 2A)(^ + Aa 



(3) 



(Bl) 
(B2) 



(B3) 
(B4) 



(B5) 



(B6) 
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